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ABSTRACT 

We report the discovery and detailed monitoring of X-ray emission associated 
with the Type lib SN2011dh using data from the Swift and Chandra satellites, 
placing it among the best studied X-ray supernovae to date. We further present 
millimeter and radio data obtained with the SMA, CARMA, and EVLA during 
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the first three weeks after explosion. Combining these observations with early 
optical photometry, we show that the panchromatic dataset is well-described by 
non-thermal synchrotron emission (radio/mm) with inverse Compton scattering 
(X-ray) of a thermal population of optical photons. We derive the properties of 
the shockwave and the circumstcllar environment and find a time-averaged shock 
velocity of v ~ 0.1c and a progenitor mass loss rate of M « 6 x 10 -5 M 0 yr _1 
(wind velocity, v w = 1000 km s^ 1 ). We show that these properties are consistent 
with the sub-class of Type lib supernovae characterized by compact progenitors 
(Type cllb) and dissimilar from those with extended progenitors (Type ellb). 
Furthermore, we consider the early optical emission in the context of a cooling 
envelope model to estimate a progenitor radius of R * « 10 11 cm, in line with the 
expectations for a Type cllb supernova. Together, these diagnostics suggest that 
the putative yellow supergiant progenitor star identified in archival HST obser- 
vations is instead a binary companion or unrelated to the supernova. Finally, we 
searched for the high energy shock breakout pulse using X-ray and gamma-ray 
observations obtained during the purported explosion date range. Based on the 
compact radius of the progenitor, we estimate that the shock breakout pulse was 
detectable with current instruments but likely missed due to their limited tempo- 
ral/spatial coverage. Future all-sky missions will regularly detect shock breakout 
emission from compact SN progenitors enabling prompt follow-up observations 
of the shockwave with the EVLA and ALMA. 

Subject headings: supernovae: specific (SN 2011dh) 


1. Introduction 


A key open question in the study of core-collapse supernovae (SNe) is the nature and 
diversity of their progenitor systems. High-resolution optical imaging of nearby galaxies 


has firmly established that the progenitors of Type IIP SNe are red super gi ants (jSmartt 


2009h. while pre-discovery imaging for T ype II SNe 1987A (Gilmozzi et al. 1987), 1993 J 
(lAldering. Humphreys fe Richmondlh 99411 and 2005gl flGal-Yam et al. 1120071: lGa.l-Ya.rn &; Leonard 
20091 ) point to a diverse set of massive stars. An alternative route to the nature of the progen- 
itors is to obtain panchromatic follow-up observations within days of explosion. The shock 
breakout pulse and subsequent adiabatic cooling of the ejecta can yield information on the 


progenitor si ze since the c 


uration and energy o 


genitor star ( 

Colgate 

1974; 

Ensman & Burrows 

1992: A 

Yaxrnan. Meszaros & Canmana 

2007; 

Chevalier & Fransson 

2008 

Katz. Budnik & Waxman 

2010; 

Nakar & Sari 

2010). Comple- 


these signals scale with the _sizej)f the pro 
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mentary follow-up observations at radio, millimeter, and X-ray bands provide unique di- 
agnostics on the shockwave velocity which scale inversely with the progenitor radius. The 
utility of such a multi- wavelength technique was demonstr ated by th e ser endipitous X-ray 


discovery and comprehensive follow-up study of SN 2008D (ISoderberg et al. 2008) 


On May 31.893 UT amateu r astro nomer Amadee Riou di scove red an optical transient 


in M51 (d ps 8.4 ± 0.6 Mpc; Feldmeie r. Ciardullo & Jacob y 1997 ). Multiple individuals 


(Griga et al. 

2011 

; Silverman, Filippenko & Cenko 

2011) 

(PTF; 

Law et al. 

2009|) reported a deep non-detection in j 


onset of the optical emission to be May 31.275-31.893 UT (Arcavi et al. 2011). Based on an 
initial spect rum on June 3.3 UT, th e trans ient was classified as a Type II supernova, dubbed 
SN2011dh (Silverman. Filippenko & Cenko 2011b Further spectroscopy reveal ed ev idence 


for h elium absorp tion features prompting the re-classification as Type lib flArcavi et al. 


2 01 ll: iMar ion et a l. 2011 ) . 


A putative progenitor star has been identified in pre-explosion Hubble Space Tele- 
scope (HST) images with a spectral energy distribution consistent with a yellow supergiant 
fiVan Dvk et al. 201 ll; Mai m d et al. 2011). The mass of the star is estimated to be between 
Mzams ~ 13 and 21 M 0 , but temperature-dependent bolometric corrections, discrepancies 
b etwee n evolutionary tracks, and treatments of rotation should also be carefully considered 
flDrout et al. 1120091) . Based on the estimated luminosity and te mperature of the object , the 


stellar radius i s R* ~ 10 13 cm flPrieto fe Hornochll201ll) . Both lYan Dvk et al. (120 llh and 


Maun d et al. (2 011 ) discuss the possibility that the yellow supergiant is instead the binary 
companion to the SN 2011dh progenitor star. In this scenario, the actual progenitor star 
may be of smaller radius. 


Recently, IChevalier fe Soderbera 020 101 ) proposed that SNe lib may be divided into two 
sub-classes based on the radius and mass loss history of the progenitor star and the properties 
of the shockwave. In this framework, compact progenitors (i?* ~ 10 11 cm) with modulated 
radio light-curves emission and shock wave v elocities of v ~ 0.1c are identified as SNe cllb 


with membe rs including SNe 200 1ig (IRvder et al. II2004I ). 2003bg (ISoderberg et al. 1120061) . 


and 2008ax ( IR.oming et al. 2009). Meanwhile, extended progenitors (77* ~ 10 13 cm) with 


smooth radio light - curves and slower s hockwaves are identified as SNe ellb (e.g., SN 1993J; 


Bartel et al. 120021; IWeiler et al. 1 120071) . The significance of the modulated radio emission 


points to an unusual (perhaps episodic) mass loss that may be unique to SNe cllb. 

SN 2011dh showed an initial peak magnitude of M n ps —16.5 mag a t At ~ 1 day since 


explosion before fading quickly flArcavi et al. 2 01 ll: [Prieto et al. 2011b The SN then re- 
brightened at At « 5 days. The two light-curve components may be interpreted as cooling 
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envelope emission followed by the radioactive decay of 56 Ni. Based on a simple compari- 
son of the SN 2011dh light-curve with that of SN 1993J and a measurement photospheric 
temperature, Arcayi et al. (20111) proposed that SN 2011dh belongs to the Type cllb class. 


Here we report the discovery and monitoring of X-ray emission associated with SN 201 ldh 
and present radio and mm-band detections from the first few weeks after explosion. We show 
that the radio and X-ray properties are consistent with those of SNe cllb and dissimilar from 
those of SNe ellb indicating that the progenitor was compact at the time of explosion. This 
is supported by our modeling of the early cooling envelope emission, which points to a pro- 
genitor radius of A* ~ 10 11 cm. Together, these diagnostics suggest that the putative yellow 
supergiant progenitor is instead a binary companion or unrelated to the SN. Finally we 
present a detailed compilation of X-ray and gamma-ray observations from multiple satellites 
and instruments obtained during the purported explosion date range. We estimate that the 
shock breakout pulse was detectable with current high energy instruments but likely missed 
due to their limited temporal/spatial coverage. 


2. Observations 

Following the optical discovery of SN 201 ldh, we initiated a prompt panchromatic follow- 
up campaign to map the non-thermal properties of the ejecta. 


2.1. Swift/XRT Observations 


Through our Target-of- Opportunity request, Swift/XRT began observing SN 201 ldh on 
une 3.50, just ~ 3 days after the explosion. As initially reported in Margutti fe Soderberg 
(2011), we discovered a bright X-ray source (S/N~ 10) at coordinates, RA J2 ooo = 13 h 30 m 5.18 s , 


Decj 2 ooo = +47 o 10'11.14" (uncertainty 4" radius, 90% confidence) at 1" from the optical SN 
position with a count-rate of ~ 0.015 cps. We analyzed 69 ks of archival pre-SN Swift/XRT 
observations and these data reveal no X-ray source at the SN position with a 3 a upper limit 
of 5.6 x 10 4 cps. This fact, coupled to the spatial coincidence of the SN, strongly suggests 
that the new source represents the X-ray counterpart to SN 2011dh (Figured]). 


Observations of SN 2011dh with Swift continued for the next several weeks. We retrieved 
and analyzed the XRT data from the HEASARC archive collected in the time period, June 
3 - July 3 UT (total exposure time of 137 ks). All XRT data were analyzed with the hea- 
SOFT (version 6.10) software package and corresponding calibration hies; standard filtering 
and screening criteria were applied. Due to the proximity of a nearby, steady X-ray source 
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(Figure [T]), we adopted a 12-pixel (~ 28") extraction region centered on the optical posi- 
tion; for lower count rates (< 0.0025 cps) we reduced the extraction region to a radius of 
6-pixels to increase the S/N ratio and eliminate contamination from a nearby faint source. 
The background was estimated from the pre-explosion Swift/X RT data to properly account 
for the contamination from the extended X-ray emission associated with M51. A spectrum 
extracted over June 3-17 UT can be modeled by an absorbed power-law with photon index 
T _ 1.5T 0.2 (90% c.l.) assuming a Galactic foreground column density N H = 1.81 x 10 20 


cm 


-2 


( Kalberla et al. 2 0 05) and no intrinsic absorption (y 2 /dof = 70.7/73, P-val=0.56). Alter- 
natively, a thermal plasma spectral model with best-fitting temperature kT = 7. 5^33 keV 
(90% c.l.) can adequately represent the data (y 2 /dof = 68.8/73, P-val=0.62). Both models 
give an average unabsorbed flux of F\ ~ 1.65 x 10 -13 erg cm -2 s -1 (0.3-8 keV) corresponding 
to a luminosity, Lx ~ 2 x 10 39 ergs _1 . 


Our resulting X-ray light-curve is shown in Figure El Over the first 10 days the SN 
faded by a factor 8. Adopting a simple power-law model for the decay, we derive an index 
of a = —0.8 ± 0.2 (90% c.l.). The source shows some evidence for spectral softening with 
time, with the photon index evolving from Tx = 0.9 ± 0.3 (90% c.l., At = 3 — 7 days) to 
r 2 = 1.8 ± 0.2 (90% c.l., At = 7 — 17 days). In comparison with X-ray observations of 
SN 1993J the softening observed for SN 2011dh begins at an earlier epoch. 


2.2. Chandra Observations 


We supplement our Swift/X RT light-curve with two Chandra observations. As reported 
byl Poolevet al. 201 ll. SN 2011dh was observed with Chandra for 10 ksec beginning on June 
12.3 UT. The SN was detected with a count rate of ~ 0.012 cps. Adopting a power-law 
spectral model, they derive an X-ray flux of F x = (1.0 ± 0.3) x 10 13 erg cm 2 s 1 (0.5-8 
keV). 


On July 3.4 UT, we obtained a second 10 ksec observation of SN 2011dh with the Chan- 
dra Advanced CCD Imaging Spectrometer (ACIS) under a Target-of- Opportunity program 
(PI Soderberg). Data were reduced with the CIAO software package (version 4.3), with cali- 
bration database CALDB (version 4.4.2). We applied standard filtering using CIAO threads 
for ACIS data. We clearly detect a source at the SN position. Extracting within a 15 pixel 
aperture, we derive a (background-subtracted) count rate of ~ 0.0049 cps. Adopting the 
same power-law spectral model described in 1 12.11 we derive an unabsorbed X-ray flux of 
F x = (2.75 ± 0.55) x 10~ 14 erg s cm 2 (0.3-8 keV). We further note that we do not detect 
any bright nearby sources that would otherwise contaminate the extraction region adopted 
for the Swift/X RT data. 
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In Table 1 and Figure [2] we report the combined Swift/XBT and Chandra light-curve 
of SN 2011dh, representing the best-sampled X-ray light-curve for a SN lib to date. We 



Figure, the SN 2011dh X-ray light-curve is more closely related to Type cllb explosions and 
a factor of ~ 10 less luminous than that observed for the Type ellb SNe 1993J and 2001gd. 
We therefore suggest that the Type cllb sub-class may be further characterized by low X-ray 
luminosities emission and early spectral softening. 


2.3. CARMA Observations 


We observed SN 2011dh with the Combined Array for Research in Millimeter- wave As- 
tronomy (CARMA; Bock et al 120061) beginning on 2011 June 4.1 UT. Observations were 
conducted with CARMA’s nine 6.1-m antennas and six 10. 4- m antennas in the D configura- 
tion, with a maximum baseline leng th of 150 m. We implemented radio and optical pointing 
flCorder. Wright fc Careen Aoiol i. We selected central frequencies of v = 100 GHz and 230 
GHz, with a total bandwidth of ~ 6 GHz. Gain calibration was performed with J1153+495 
and we used a source-calibrator cycle time of ~ 15-20 minutes. Flux and bandpass calibration 
was carried out using observations of 3C273 and 3C345 and Neptune resulting in an overall 
uncertainty in the absolute flux calibration of ~ 10%. We used t he Multichannel Image Re- 
construction Image Analysis and Display (MIRIAD; Saul t. Teuben & Wright 19950 software 
package for data reduction. We integrated on SN2011dh for 43 minutes at each frequency 
and clearly detected a source coincident with the optical and X-ray SN positions. Preliminary 
results for the v = 100 GHz observation were presented by IH oresh . Zauderer & Car p enter 
(2011). Here we present the results of the v = 230 GHz observation and a re-analysis of the 
v = 100 GHz data in which we have refined the flux calibration. Fitting a Gaussian model 
to the source, we derive an integrated flux density of F u = 4.5 ± 0.3 mJy at v = 100 GHz 
and a 3cr upper limit of F u < 3.5 mJy at v = 230 GHz (Table 2). Thus the radio spectrum 
is optically-thin between the two mm-bands at At « 4 days since explosion. 


2.4. SMA Observations 


Contemporaneously with the CARMA o bserv ations, we observed SN 2011dh with the 
Submillimeter Array (SMA; Ho, Moran & Lo 2004) on June 4.0 UT in the compact configu- 
ration at a frequency of v = 230 GHz with 8 GHz bandwidth. Observations included all eight 


antennas. Passband calibration was performed in the standard way using Neptune, 3C454.3 
and 3C279. We used 3C279 for flux calibration, verifying our calibration with observations 
of Titan in addition to 3C279 and our gain calibrators on June 10 UT. The absolute flux 
calibration is accurate to ~10%. We flagged low elevation data (< 21°) and the first several 
hours of the observation when weather conditions were less favorable, prior to improved 
pointing solutions, and when one antenna was missing. The data were calibrated using stan- 
dard MIR/IDL routines developed for the SMA, with further calibration and im aging c arried 
out in MIRIAD and the Astronomical Image Processing System (AIPS; Greisen 120031) . The 
resulting total integration time on source was 2.75 hours. We detect a radio source coincident 
with the SN position with a flux density of F u = 3.6 ± 0.9 mJy (Table 2). This detection is 
consistent with the 3cr upper limit from CARMA. 


2.5. EVLA Observations 


Array (EVLA: 

Perlev et al. 

2009 

with a flux density of F u ps 2.68 ± 0.10 mJy at v = 22.5 

GHz f 

Horesli et al. 

2011 

). There 

is no coincident radio source in the catalog of M51 compact 

radio 

sources ( 

Maddox et al. 

2007 

). A comparison with the CARMA and SMA flux densities 


obtained contemporaneously indicates that the spectral peak lies between the EVLA and 
CARMA bands (Figure E|). 


We began monitoring SN2011dh with the EVLA on June 17 UT (At m 17 days after 
explosion) as part of a Rapid Response observing program for long-term monitoring of the 
supernova (PI Soderberg). Data were collected within the wide C, X, Ku, K, and Ka bands. 
Within each of these bands (except X) we selected two central frequencies enabling spectral 
coverage spanning v = 5.0 — 36.0 GHz. Each central frequency has associated bandwidth of 
0.8 to 1.0 GHz. All observations were obtained in the (most extended) A-array. We used 
J1327+4326 to monitor the phase while the absolute flux calibration was carried out using 
3C286. Data were reduced using AIPS and the Common Astronomy Software Applications 
(CASA). We fit a Gaussian model to the radio SN emission in each observation to derive 
the integrated flux density (Table 2). The reported flux density errors include errors from 
Gaussian fitting, the map rms noise, and systematic errors of 1% at low frequencies (5-16 
GHz) and 3% at high frequencies (20.5-36.0 GHz). At this epoch, the peak of the radio 
spectrum has clearly shifted to the cm-band. 


Additional observations with the EVLA are on-going and are the focus of a separate 
paper. In Figure [3] we compare the radio spectrum of SN 2011dh with the newly available 
and unprecedented wide bands of EVLA which enable continuous spectral coverage from 
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~ 1 — 40 GHz. SN2011dh represents the first SN for which such detailed mapping of the 
spectrum has been possible in the EVLA era. 


A Model for the Radio Emission 


Early radio observations of SNe uniquely trace the shockwave as it races ah ead o f the 
bulk ejecta and shock- accelerates particles in the local circumstellar medium (CSM; Chevalier 
19821 ). This environment was enriched by the progenitor star wind during the centuries 
leading up to the explosion. Through this dynamical interaction, the shockwave accelerates 
CSM electrons into a power- law distribution, A(y) oc r )~ p , above a minimum Lorentz factor, 
7 m . The accelerated electrons gyrate in amplified magnetic fields giving rise to non-thermal 
synchrotron emission. In the case of Type Ibc and cllb supernovae, the radio emission is 
quenched at low frequencies primarily due to synchrotron self-absorption (SSA), producing 
a spectral turnover that defines the peak of the radio spectrum, v p . The self-absorbed radio 
spectrum is described by F v oc v h ^ 2 below v p and F v oc i/Mp- 1 )/ 2 above u p . As shown i 


m 


Figure [3j our EVLA, CARMA, and SMA observations of SN 2011dh on two separate epochs 
are well-described by a synchrotron self-absorbed spectrum with p « 3. We note that the 
modest disagreement between the SSA model and the measurements near the spectral peak 
may indicate asphericity of the emitting region. 


Chevalier (119981 ) showed that for radio SNe with minimal free-free absorption, the radius 
of the shockwave, R, and its time-averaged velocity, v, can be robustly estimated from the 
observed values of u p and the associated peak spectral luminosity, L UtP . For p « 3, the shock- 
wave radius is given by R ~ 3.3xl0 15 (e e /e j B) _1 / 19 (L !/p ? «) 8 / 1 v r i cm where L, y „ ^ i s normalized 


to 10 26 erg s 1 Hz 1 and z/ Pj5 is normalized to 5 GHz dChevalier fe Franssonll2006l ). The frac- 


tions of post-shock energy density shared by accelerated electrons and amplified magnetic 
fields are denoted by e e and e#, respectively, and we assume that the radio emitting region is 
half of the total volume enclosed by a spherical shockwave (e.g.. Chevalier & Fransson 2006). 


As shown in Figure [21 for SN 2011dh we find v p ps 40 and ~ 11 GHz on 2011 June 4 and 
17 (A t ~ 4 and 17 days), respectively, with an associated peak luminosity of and L v 2 q ~ 6.7 
on both epochs. These observables correspond to shockwave radii of R ~ (1.0, 3.7) x 10 15 
cm assuming typical partition values of e e = es = 0.1. The time averaged shockwave velocity 
is thus v ~ 0.1 c and thus a factor of ~ 2 faster than the material at optical photosphere 
at At ps 3 days (Silverman. Filippenko & Cenkol 2011 ). The total internal energy required 
to power the observed radio signal can be est imated from th e pos t-shock magnetic energy 
density, E = B 2 R^ /12eB- As shown by Chevalier & Fransson (120061 ). the amplified magnetic 
field is directly determined from the spectral properties, B pa 0.70 (e e / eg) -4 / 19 L^g 9 v p $ G. 
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At At m 4 and 17 days, we find B m 4.5 and 1.2 G. The total internal energy is thus, 
E m (1.7, 6.3) x 10 46 erg for the two epochs, respectively, by maintaining the assumption 
that €b = 0.1, i.e. 6b, - i- The roughly linear temporal increase of internal energy is consistent 
with a slightly decelerated shockwave. As shown in Figure 01 the shockwave properties of 
SN 201 ldh are similar to those of Type cllb, which tend to also be characterized by shockwave 
velocities of v > 0.1c while slower shockwaves are inferred for SNe ellb. 


-l 


The progenitor mass loss rate is M m 0.39x 10 -5 es 1 _i(e e /cs) - 8 ^ 19 -^^ 26 9 ^, 5^,10 y r 

where f Pil0 is the observed time of the spectral peak normalized to 10 days and we have 
assumed a wind velocity of v w = 1000 km s _1 . Thus, we estimate a progenitor mass loss 
rate for SN2011dh of M ~ 3 x 10 -5 M 0 yr -1 , similar to the mass loss rates derived for 
SNe Ibc and cllb and also similar to the values ob served for Galactic Wolf-Rayet stars 
(iCappa. Goss fe van der Huchtl 2004 ; Crowther 2007 ). We note that in this framework, the 
radio data constrain the ratio, (M/v w ), such that a variation in the assumed value of v w 
shifts the mass loss rate estimate by the same factor. 


4. A Model for the X-ray Emission 


On timescales of days after explosion, the X-ray emission observed from SNe Ibc and 
lib may be dominated by a number of different emission processes including synchrotron, 
thermal, and inverse Compton scattering. Continued energy input from a compact remnant 
(black hole, magnetar) has also been invoked to explain the X-ray emission from som e 
events (e.g., SN2006aj, Soderberg et al. 2006 ; SN 1979C, Patnaude. Loeb fe Jones 2011 b 
The temporal and spectral evolution of the radio and X-ray emission, together with their 
observed luminosity values, enables these different processes to be distinguished. Here we 
consider the nature of the observed X-ray light-curve for SN 2011dh. 

The radio-to- X-ray spectral index is observed to be (3rx ~ —0.7 on both June 4 and 
17 UT. Thus, an extrapolation of the optically-thin radio synchrotron spectrum with pm 3 
under-estimates the X-ray flux at both epochs by a factor of ~ 140. Sim ilarly high radio- 


to- X-ray spectral indices are not atypical for SNe Ibc and lib (see Cheval i er fe Franss on 
2006 and references within). Attributing the X-rays to synchrotron emission would require 
efficient particle acceleration and a flattening of the electron distribution at high energies. 


As shown in I Cheval ie r fe Frans son 020031 ). the thermal emission from shock heated ejecta 
in the reverse shock region may give rise to strong free-free X-ray emission. If the material is 
fully ionized and in temperature equilibrium, the expected luminosity is broadly determined 
by the electron temperature and the mass loss rate of the progenitor star with an expected 
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linear decay in time. We use equations (25-30) of Chevalier & Fransson (2006) with v ~ 0.1c 
and M _5 fa 3. We estimate the free-free X-ray luminosity for SN2011dh to be Lff & 
2.5 x 10 37 erg s _1 on June 4 UT which is a factor of ~ 200 lower than the Swift/XRT 
measurement at this epoch (Figure [2]). 

We next consider an inverse Compton scattering model in which the optical photons 
associated with both envelope cooling emission and the radi oac tive decay of 5 6 Ni are upscat- 
tered to the X-ray band by radio emitting electrons (j Chevalier. Fransson & Nvmark 2006 ). 
In this scenario, the X-ray decay should track the optical evolution through the cooling en- 
velope decay to the re-brightening due to 56 Ni decay. The early optical emission indicates 
an average luminosity of L boi ~ few x 10 42 erg s ^ 1 and a minimum near At & 5 days. The 
X-ray light-curve suggests a similar minimum near this epoch and an overall decay of roughly 
Lx oc t~ l (Figure [ 21 ) • Adopting the formalism of [Chevalier & Fransson! (12006! ) and assuming 
the relativistic electron population extends down to (y m = 1), the predicted IC emission is 
Lie ~ 2.9 x 10 36 es -iM-5(e e /eB) n / 19 (u/0.1c) -1 Lboi,42At^ 0 erg s _1 where M _ 5 is the mass 
loss rate normalized to 10 ~ 5 M 0 yr ” 1 and we maintain the assumption of v w = 10 3 km s _1 . 
Here, 10 is the time since explosion normalized to 10 days. 

We find that the X-ray emission may be attributed to inverse Compton emission if the 
assumption of equipartition is relaxed and (e e /es) ~ 30 with es ~ 0.01. This deviation from 
equipartition implies modest adjustments to our physical parameters estimates, including a 
factor of ~ 2 increase in the mass loss rate and a factor of ~ 4 increase in the internal energy 
of the radio emitting material. The modified values for the mass loss rate and magnetic held 
are M _ 5 rs 6 and 5 r 0.4 G, respectively. A prediction of this model is that the X-ray light- 
curve will decay more steeply following the optical SN peak ([Chevali er & Fransson 20061). 
which may be suggested by our Chandra measurement on July 3 UT. 


5. Constraints on the Progenitor Size 


The early optical emission from SNe is dominated by the adiabatic cooling of envel ope 
mate rial following the breakout of the shockwave through the stellar surface (lEnsman fe Burrows 
1992 ). This component cascades through the optical band in the first few days following 
explosion. The radius and temperatur e assoc i ated with this component may be roughly 
approximated as thermal (although see Nakar & Sari 2010 for a more comprehensive discus- 
sion of the spectral evolution) and used together with estimates for the bulk SN parameters 
(including the ejecta kinetic energy, Ek, and mass, M e j) leads to a determination of the pro- 
genitor radius, A*. Such techniques have been used to derive the progenitor radius of SNe 
1987A, 1999ex, 2008D, and 2008ax in addition to SN 2006aj associated with XRF 060218. 
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Here we adopt the formalism of Chevalier & Fransson (2008J) in which the temperature of 
the photosphere is given by T ph ft 7800 E [) K { S X Mf^f/ 4 R ( f 2 x /Atf 0AS K and the photospheric 
radius is R p h ~3x 10 6 14 * E/f^ M e j° 0 28 * At°- 78 cm. Here we have normalized E K5 i to 10 51 erg, 
M e j i0 to M 0 and R*,n to 10 11 cm. Within this framework, it is assu med that the density 


and pressure profiles of the envelope are consistent with those of Matzncr & McKee (119991) . 
and that the bolometric luminosity decays as L ph oc A tf ou . 

We model the first few optical obse rvat ions (At A 5 da ys) including the i nitial detect ion 
of Lph 10 42 * erg s _1 on June 1.191 UT ( Arcavi et al. 1 2 0111 ). Silv erma n. Filippenko & C enkc 
(2011) report a photospheric velocity of v p h ~ 17, 600 km s _1 at At ft 3 days after explosion 
implying a ratio of (Ek, 51 / M^®) ~ few. We estimate R p h ft 5 x 10 14 cm and T p h ft 8000 K at 
At ft 1 day. This is rough ly con sistent with the photospheric temperature derived from op- 
tical spectroscopy flArcavi et al. 20111). and implies a compact progenitor since T ph oc R °' 25 
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The high luminosity of the initial detection requires a ratio, Ek/M s j ~few, in line with the 
ratio implied by the high photospheric velocity. Thus, the early optical emission points to 
a compact progeni tor, similar to those of SNe Ibc and cllb, and consistent with the earlier 
report by A rcav i et al. (2011). 

We note, however, that the observed optical decay, L oc A tf 1 , is significantly steeper 
than the model prediction. This may point to an irregular density profile near the stellar 
surface, perhaps associated with an unusual mass loss ejection in the final stage of the 
progenitor’s evolution. Such an irregular density profile may affect some of the conclusions 
that we derived above based on the canonical model and this will be the focus of a future 
publication. We note that unusual mass loss histories have similarly been inferred for other 
SNe cllb based on radio observations (e.g., SN 2003bg; ISoderberg et al. 2006 ). 


6. Shock Breakout X-ray Emission 


For compact progenitors such as SN 2011dh, the b reako ut pulse may be detectable at 


X-ray and gamma-ray energies (e.g., SN 2008D; ISoderberg et al. 120081 ). To this end, we 


searched for evidence of a high energy pulse associated with the shock breakout of SN 2011dh 
using data collected by the Swift Burst Alert Telescope (BAT), the Monitor of All-sky X-ray 

Image (MAXI) camera attached to the Japanese Experiment Module, and the Interplanetary 

Network (IPN: Mars Odyssey, Konus-Wind, RHESSI, INTEGRAL (SPI-ACS), Swift /BAT, 

Suzaku, AGILE, and Fermi/GBM; we note that MESSENGER was in superior conjunction 

and off during this period). 


The position of SN2011dh was observed by Swift / BAT over many pointings during 
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the estimated explosion date range, May 31.275-31.893 UT. No gamma-ray emission was 
detected from the SN. We compile the Swift/ BAT observations in Figure 0 Each pointing 
had a duration of ~ 500 — 1000 sec, and the total time on source was 11148 sec or 20% of 
the full explosion date range. We estimate the count rate at the SN position and also in the 
background in each pointing and note that the sensitivity varies as a function of the off-axis 
angle, ranging from 10 to 57 degrees. Assuming a power law spectrum for the breakout flux 
with a photon index of T = 2, we infer 5cr upper limits on the gamma-ray emission from the 
SN of F 7 < (1.1 — 3.9) x 10” 9 erg cm -2 s” 1 (14-195 keV) for each individual pointing. The 
combined upper limit is F 1 & 2.8 x 10” 10 erg cm -2 s” 1 for the same energy range. At the 
distance of M51, this limit corresponds to a luminosity of L 7 < 2.4 x 10 42 erg s” 1 . 


MAXI flMatsuoka et al. 1 120091 ) is an X-ray all sky monitor that scans most of the sky 
every 92 minutes with its Gas Slit Cameras (GSC; Mihara et al. Eon]). Unfortunately during 
the estimated period of the SN breakout, the direction of M51 was covered only by a camera 
with high background, resulting in less stringent limits than usual. The X-ray image taken at 
every scan transit (~ 80 sec), encoded in one dimension by the detector position and the other 
by transit timing, was fitted to the PSF with the fixed source position to evaluate the source 
count rate. The energy range of (4-10) keV was chosen to maximize the signal-to-noise ratio. 
Then the count rate was converted to energy flux assuming a power-law model with a photon 
index of T = 2 resulting in a conversion factor of F\ ~ 1.1 x 10” 8 erg cm -2 s” 1 . We derive 
an average 5 a upper limit of Lx < 3.9 x 10 42 erg s” 1 . In comparison, the observed breakout 
pulse from SN 2008D ( d m 30 Mpc) had a luminosity ~ 25 times higher (Lx 
0.3-10 keV) with a duration of ~ 5 min. 


10 44 erg s 4 ; 


While Swift /BAT and MAXI both have a limited temporal coverage, the IPN is full 
sky with temporal duty cycle_of nearby 1 00%; it is sensitive to hard X-ray emission in the 
energy range, (25-150) keV flHurlev et al. 20 10). Within a two-day window centered on the 
estimated explosion date, a total of four triggered bursts were detected by the IPN. All four 
were confirmed through observations by multiple instruments or spacecrafts and thus could 
be localized. These detections include Fermi/GBM (3 bursts), Konus (1), INTEGRAL (1) 
and Swift / BAT (1). In all cases the localizations were statistically inconsistent at a > 3 a level 
with the position of SN 201 ldh. Thus we find no evidence for a detected gamma-ray transient 
in coincidence with SN 2011dh and adopt a hard X-ray upper limit of F < 6 x 10”' erg cm” 2 


corresponding to an energy of 7% < 5 x 10 45 erg for SN 201 ldh. 
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6.1. A Comparison to Breakout Predictions 


Here we compare theoretical predictions for a prompt breakout pulse with the avail- 
able X-ray and 7 — ray limits. We consider a model in which the progenitor is compact, 
iL n ~ 1, and emb edde d in a stellar wind with M _ 5 « 6 . In the standard SN Ibc model 
of Chevalier & Fransson (120061 ), also applicable to SNe cllb, the shockwave radius evolves 
as R oc t m with m ~ 0.9. Extrapolating back to the breakout radius, the shockwave speed 
at the stellar surface is n 0 ~ 0.3c and thus /3q ~ 0.3/?_o.s. In the case of such fast shock 
velocities, there are deviations from thermal equilibrium resulting in high energy emission 
( Katz, Budnik fc Waxman 2010 ; Nakar fe Sari 2010 : Sapir, Katz fe Waxman 2011 ). 

The optical depth at the stellar surface is r = kM /A nv w r ~ 2M_ 5 i?b 1 1 where we adopt 
k = 0.4 cm 2 g _1 and Ru is the shock radius normalized to 10 11 cm. The shock breakout 
will occur outside of the star, within the stellar wind at a radius, i?br = /3qkM / 4ttv w ~ 
6 x 10 10 M_ 5 /3_o.5 cm as long as R * < R^. I 11 the case of SN 2011dh, we therefore predict 
a wind breakout at radius, R}, r ~ 4 x 1 0 11 cm. The_energy associated with breakout is 
Eb r k 3 x 10 43 M 2 5 /9 3 0 5 ( Katz. Sapir fe Waxman 2011). On the breakout timescale, a 
radiative collisionlcss shock forms and an additional comparable or larger amount of energy is 
emitted (i Kat z. Sapir & Waxman 2011). Adopting our parameters for SN 2011clh, we predict 
Ebr > 10 45 erg. During the transition from a radiation mediated shock to collisionless, the 
temperature rises steadily from > keV to > 100 keV and we roughly estimate the spectrum 
to be uF u ~ constant. 


The predicted breakout energy is within the detectability range of the Swift/BAT and 
MAXI upper limits spanning (4-195) keV, however, both offer only limited temporal coverage. 
We estimate the rise time of the breakout pulse to be tb r ~ Rb r /v 0 & 1 min. Efficient 
emission from the collisionless shock continues beyond this time, and the full duration of the 
pulse is necessarily longer. Given the large gaps in the Swift / BAT and MAXI coverage (see 
Figure [5]), especially toward the beginning of the explosion date range, we speculate that 
Swift and MAXI missed the breakout pulse. The IPN upper limit is a factor of a few above 
the predicted energy and thus consistent with this breakout model. Finally we note that 
if the progenitor was more extended with a radius, R * ~ 10 13 R 13 cm, a stronger pulse of 
Ebr ~ 3 x 10 4 ' Rf 13 erg is expected at lower frequencies, hv < 300eV. Such a pulse is not 
constrained by the hard X-ray/gamma-ray observations presented here. 
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7. Conclusions 

We present multi- wavelength follow-up observations of SN2011dh spanning the radio, 
millimeter, X-ray, and gamma-ray bands and obtained within the first few seconds to weeks 
following the explosion. The X-ray light-curve for SN 2011dh is perhaps the best-sampled 
to date for a SN lib and suggests that X-ray properties (luminosity, spectral evolution) may 
further distinguish compact and extended progenitors. Using the newly available wide bands 
of the EVLA, we show that the radio and millimeter data are consistent with a synchrotron 
self-absorbed spectrum, while the X-ray emission may be understood as inverse Compton 
upscattering of optical SN photons. These data offer unique diagnostics (shockwave velocity, 
v ~ 0.1c) and point to a compact progenitor star. Through modeling of the early optical 
emission we find that the photosphere is characterized by a low temperature, T p h m 8000 K, 
and the physical parameters of the ejecta are constrained to be E K 51 /M e ^ Q ~ few. These 
properties point to a compact progenitor, R * ~ 10 11 cm, however, we find that the fast decay 
of the early optical emission at At < 5 day is incompatible with canonical cooling envelope 
models and may suggest an irregular ejecta profile. A compact progenitor size is inconsistent 
with the extended radius (/?,* ~ 10 13 cm) of the coincident yellow supergiant identified in 
pre-explosion HST imaging; we suggest that this object is either a binary companion or 
unrelated to the SN. We conclude that SN 2011dh is a likely a member of the Type cllb class 
of core-collapse explosions. Our long-term EVLA monitoring observations will reveal if the 
radio light-curves are modulated, indicative of a variable and/or episodic mass loss history 
in the decades leading up to explosion - an observational characteristic shared by other SNe 
cllb. 


Finally, we used the parameters derived above for the shockwave and progenitor to 
estimate that the shock breakout pulse from SN2011clh was detectable in the X-ray and 
gamma-ray bands with current satellites. We attribute the lack of a detection as likely 
due to limited temporal/spatial coverage during the estimated explosion date. Looking 
forward, sensitive and wide- field X-ray experiments such as LOBSTER and JANUS will 
regularly di scove r sh ock b reakout emission from similarly nearby and compact SN progenitors 
Soderberg et al. 2009 for a discussion). Such detections will not only provide 


tsee e.g. 

information on the progenitor size, but also accurate explosion date estimates will enable 


statistically significant search es for gravitational wave and neutrino counterparts (IOttll20Q9 


Katz, Sapir & Waxman 2011)- Furthermore, these prompt discoveries will enable rapid low- 


frequency follow-up with sensitive radio and millimeter facilities thanks to the advent of 
EVLA and the Atacama Large Millimeter Array (ALMA). 
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Fig. 1. — (a): A (0.3-10) keV co-added archival image (total exposure of 64 ks) of M51 
obtained pre-explosion using Swift /X RT in the time range June 2006 - May 2009. No sta- 
tistically significant excess is detected at the SN2011dh position, (b): XRT observations 
of M51 (total exposure of 78 ks) obtained after the discovery of SN 2011dh clearly reveal 
an X-ray source at the SN position, (c): Archival Chandra observations reveal no source 
within 28 arcsec of the SN position (circle), (d): SN2011dh is detected in our Chandra 
ToO observation on July 3. No bright contaminating X-ray sources are detected within the 
extraction region. We restrict our XRT extraction region to 6-pix at late times to avoid the 
faint source to the SW. 
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Fig. 2. — X-ray emission from SNe 

lib: SN1993J 

(Chandra et al. 

2009); SN 2001gd 

(Perez-Torres et al. 

2005); SN 2008ax ( 

Roming et al. 

2009 

); SN 2001ig 

(Schlcgcl & Ryder 

2002 

), and SN 2003bg ( 

Soderbere et al. 

2006) and SN 201 ldh (this work) 

Error bars are ler 


For SN 2011dh the energy band (0.3-8) keV is used to allow a direct comparison to SN 1993J. 
SNe of Type cllb are shown in red while Type ellb are shown in blue. The X-ray luminosities 
of Type cllb appear lower than those of Type ellb. 
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Fig. 3. — The radio spectrum of SN 2011dh across multiple epochs - At ~ 4 (red) and 17 
(blue) days - is well described by a synchrotron self-absorbed spectral model with F oc 
(F u oc z/“( p_1 b 2 ) above (below) the spectral peak, u p . The observations indicate an electron 
energy index of p ~ 3. Error bars are la. 


23 



Fig. 4. — The radio properties of various SNe are compared including SNe lb (green) Ic 
(cyan), cllb (orange), ellb (red), and engine-driven SNe associated with nearby GRBs (blue). 
The properties of the spectral peak can be used to derive the shockwave velocity (dashed 
lines). With a blastwave velocity of v ~ 0.1c, SN 2011dh (yello w star) is more similar to 
Types Ibc and cllb than Type ellb. Adapted from [C heval ier fe So derber g (2010). 
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Date (May 31 UT) 

Fig. 5. — Swift/BAT and MAXI observations of SN2011dh during the estimated explosion 
date range, May 31.275-31.893 UT in the energy range (14-195) keV and (4-10) keV, respec- 
tively. No source is detected at the position of the SN during the individual (~ 1000 sec) 
BAT pointings or (~ 80 sec) MAXI scans. The time ranges when the SN was in the FOV of 
the instruments are shown in grey. 
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Table 1. X-ray Observations of SN2011dh 


Date 

Time Range 

Unabsorbed Flux 1 

Error 

Satellite 

(UT) 

(days) 

(erg cm -2 s^ 1 ) 

(erg cm -2 s _1 ) 


June 3.512 

0.004 

1.12 x 10- 12 

2.8 x 10- 13 

Swift/X RT 

June 3.526 

0.003 

6.28 x 10" 13 

2.8 x 10- 13 

Swift/X RT 

June 3.585 

0.007 

7.45 x 10- 13 

1.7 x 10- 13 

Swift/X RT 

June 3.943 

0.152 

5.01 x 10- 13 

2.2 x 10" 13 

Swift/X RT 

June 4.753 

0.252 

4.25 x 10" 13 

1.8 x 10- 13 

Swift/X RT 

June 5.326 

0.033 

3.80 x 10” 13 

1.1 x 10" 13 

Swift/X RT 

June 6.536 

0.571 

1.29 x 10- 13 

4.9 x 10" 14 

Swift/X RT 

June 8.042 

0.181 

1.34 x 10~ 13 

6.0 x 10" 14 

Swift/X RT 

June 8.911 

0.253 

3.69 x 10" 13 

1.1 x 10- 13 

Swift/X RT 

June 10.020 

0.301 

1.96 x 10~ 13 

6.8 x 10" 14 

Swift/X RT 

June 10.954 

0.165 

2.16 x 10- 13 

8.9 x 10" 14 

Swift/X RT 

June 12.027 

0.370 

1.87 x 10~ 13 

6.7 x 10 -14 

Swift/XKY 

June 13.127 

0.179 

2.03 x 10” 13 

7.4 x 10 -14 

Swift/X RT 

June 13.887 

0.200 

1.89 x 10- 13 

6.6 x 10 -14 

Swift/X RT 

June 15.493 

0.601 

1.80 x 10” 13 

6.4 x 10~ 14 

Swift/X RT 

June 16.969 

0.136 

1.23 x 10- 13 

4.8 x 10 -14 

Swift/X RT 

June 21.345 

2.848 

1.06 x 10" 13 

2.5 x 10" 14 

Swift/X RT 

June 29.278 

4.227 

6.57 x 10" 14 

2.2 x 10" 14 

Swift/X RT 

June 13.300 

0.115 

1.00 x 10" 13 

3.0 x 10- 13 

Chandra/ ACIS* 

July 3.500 

0.115 

2.77 x 10~ 14 

5.5 x 10" 15 

Chandra/ ACIS 


1 Energy range, (0.3-8) keV. 


t Observation reported by Pooleyef al. (2 0111 ) within energy range (0.5-8) keV. 


Tabic 2. 


Radio and Millimeter Observations of SN 2011dh 


Date 

Central Frequency 

Flux Density 

Error 

Telescope 

(UT) 

(GHz) 

(mJy) 

(mJy) 


June 4 

22.5 

2.68 

0.10 

EVLA* 


100 

4.5 

0.3 

CARMA 


230 

<3.5 


CARMA 


230 

3.6 

0.9 

SMA 

June 17 

5.0 

2.430 

0.037 

EVLA 


6.8 

4.090 

0.048 

EVLA 


8.4 

5.535 

0.057 

EVLA 


13.5 

6.805 

0.072 

EVLA 


16.0 

6.721 

0.070 

EVLA 


20.5 

6.472 

0.195 

EVLA 


25.0 

5.127 

0.155 

EVLA 


29.0 

4.603 

0.140 

EVLA 


36.0 

3.473 

0.108 

EVLA 


t Observation reported by Horesh et al 12 01111 . 


